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Abstract 

By applying fractional calculus to the equation proposed by M. Planck in 1900, we obtain a new blackbody 
radiation law described by a Mittag-Leffler (ML) function. We have analyzed NASA COBE data by means of a 
non-extensive formula with a parameter (q - 1), a formula proposed by Ertik et al. with a fractional parameter 
(a - 1), and our new formula including a parameter (p - 1), as well as the Bose-Einstein distribution with a 
dimensionless chemical potential p. It can be said that one role of the fractional parameter (p - 1) is almost the 
same as that of chemical potential (p) as well as that of the parameter (q - 1) in the non-extensive approach. 

Keywords: Planck distribution; fractional calculus; Mittag-Leffler function; Bose-Einstein distribution; 
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1. Introduction 


The NASA COBE Collaboration has reported that the universe is full of photons at a temperature of 2.725 K 
JUSSI . Their distribution is described by the Planck blackbody radiation law as follows: 

f/pianckCT, V) = 2X (1) 

e x - 1 

where Cb = 8 nhv^/c 3 , x = hv/ksT, h, ks, and T are the Planck’s constant, the Boltzmann’s constant, and temper¬ 
ature, respectively, c is the speed of the light. 

Moreover, the following residual spectrum has been reported: 

[residual spectrum] = [COBE data] - [Eq. 0} with 7 cmb K], 


It is worthwhile to notice that there are two kinds of residual spectra (1994) and (1996). To explain the resid¬ 
ual spectrum mentioned above JSSS] , the following Bose-Einstein distribution with a dimensionless chemical 
potential p is utilized by the NASA COBE Collaboration 0: 


U (BE \T , v, p) = 




e x +r - l 


-> Cb 


e x - 1 




(e x -l Y 


( 2 ) 


To investigate the residual spectrum (1994) Hi , Tsallis et al. 0] computed the following formula based on a 
non-extensive approach. 


j/NETD l \T, V, q) = £/pi anck (r, V) + C B f ~4 

e x - 1 
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ln(l - e~ x ) - x 


1 + e~ x 
1 - e~ x 


x 2 1 + 3c-* 
+ ~2 (1 - e~ x ) 2 


(3) 


where (q — 1) is known as the non-extensive parameter [00!- The correction term to the Planck distribution in 
Eq. ([3]) is carefully calculated in [0. The residual spectrum (1996) is investigated in 0. 
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Second, it is also known that Tirnakli et al. [9] have calculated the following formula based on the non-extensive 
approach, i.e., q algebra (see Ref. 0, 111 [Hill): 


^(NETDIt)^ y q) = 


Cr 


(«-!)«! 


[1 +(q- l)x]!/(«-9- 1 


■* Cr 


1 <7-1 e x x 2 

+ 


e x - 1 


2 (e x - l) 2 


(4) 


Third, using a generalized partition function derived from 114, 15, I . Ertik et al. have proposed the following 
formula based on fractional calculus with the Mittag-Leffler (ML) function El: 


U (FC l \T, v, a) = 


C B 


(or—l)<$cl 


■» Cr 


E a (x) - 1 

where the Mittag-Leffler (ML) function is defined as follows: 

x" 


1 , n /« 

+ (a- 1 )- 


e x — 1 


(e x - l) 2 


(5) 


Ax) = ^ 


n=0 


r (na + 1) 


The function f(x) is expressed as follows: 


CO 


fix) = 2 

k =0 


<bcV(l +k) 

T(1 + k) ’ 


where ip(z) = </(lnT(z))/ dz is the digamma function, and (a - 1) is known as the fractional parameter. Those 
distributions are summarized in Table Q] 


Table 1: Typical modified Planck distributions 


Bose-Einstein distributions 

I 

pX+R - 1 

/r : dimensionless chemical potential 

Non-extensive approach II 

1 

e\-\ 

q : non-extensive parameter 
e x = [1 +(1 - ? )jr] 1/(1 -«> 

Fractional calculus I 

1 

E a (x) : Mittag-Leffler function 

EJx) - 1 

a : fractional parameter 



In this study, we investigate the third approach, the fractional calculus, in more detail, because we are interested 
in the approach in El- Concerning the above-mentioned problem, we would like to adopt a different viewpoint 
from that of Ertik et al. El- It is well known that in the derivation of the blackbody radiation law in 1900, Planck 
adopted the following thermodynamical equation 1 18 . El El (See also 111): 


dU_ 

dp 


-aU -bU 2 , 


( 6 ) 


where U is the energy density, p = 1 /k B T is the inverse temperature, and a and b are parameters. Using the 
ordinary calculus, we obtain the following expression: 

U(x)= a J b , x = a(p - Pq). (7) 

e x - 1 

If we introduce a fractional derivative instead of the partial derivative for p in Eq. ©, we cannot derive 
an analytical solution, because the equation is nonlinear for the function U. Therefore, we put U — 1/R and 
x = a(P - Pa). Then, we obtain the following equation for R from Eq. ([6j: 


The solution for Eq. ® is given by 


dR b 

7T= r + ~- 

dx a 


( 8 ) 


R(x) - 
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U(xj 



(9) 
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Equation ([Sj is linear for R. Therefore, when the Riemann-Liouville fractional derivative is introduced in Eq. © 
instead of the partial derivative of x, the Mittag-Leffler function appears in the solution [21, 22, [231 24, [25]. 

In this study, we aim to apply fractional calculus to Eq. © in $2] In 5(3] various properties of the new formula 
in addition to Eqs. ([2]», ©, and © are investigated. In © analysis of the NASA COBE data in terms of the new 
formula, as well as Eqs. ©, ©- ©, and © are presented. In © the concluding remarks and discussions are 
presented. 

In |Appendix A| we explain an introduction of the fractional calculus in Eq. ©. In |Appendix B| Caputo 
derivative is mentioned. In Appendix C an interrelation between chemical potential p and fractional parameter 


(p - 1) is investigated. The same investigation in cases of NEXT II and FC I are briefly mentioned. 


2. Application of Fractional Calculus to Eq. © 

The Riemann-Liouville fractional derivative 122. 23| of function R(x) for m =1, 2, ... is defined as follows: 

1 / d \ m r x 

o D p R(x) = — -—) I (x - T)' n ~ p ~ l R(t)cIt for rn - 1 < p < m. (10) 

r(/72 — p) \ CIX / Jo 

From Eq. d), we have o D p R{x) = j-—R(x), if p = m - 1. The Riemann-Liouville fractional integral is defined 
as follows: 

0 D~ p R(x) = f (x - t) p ~'R(t)(1t for 0 < p. (11) 

r ip) Jo 

To obtain a fractional blackbody radiation formula or Bose-Einstein distribution, the partial derivative for x in 
Eq. ® is replaced by the Riemann-Liouville fractional derivative (ITOb . Then, the following equation is obtained: 

o D p R(x) = R(x) + b/a, x = a(J3 -yS 0 ) > 0. (12) 

If p = 1, Eq. ( 1 1 21 } is reduced to a first order partial differential equation, Eq. ©. Therefore, we seek the 
solution for 0 < p < 2. 

The Laplace transform R(s) of function R(x) is defined as follows: 

/-»oo 

R(s) = £[R( x y, s] = I e- sx R(x)dx. (13) 

Jo 

Using the following formulas, 


£[ 0 D p R(x); s] = s p R(s ) - J oD p - k ~ l R{x)\ x= o = s p R(s ) - £ c k s k . 


k =0 


k =0 


£[l;s] = 1/s, 


where q = oD p x k 1 /?(jc)|j C= o, we obtain from Eq. the Laplace transform R(s ), 


b s k 

R(s ) = ———— +) c k 


as(s p - 1) 4^' s p - 1 


To obtain the solution R(x) from Eq. we use the following formulas 112211 of the Laplace transform, 

L[xP~ l E afi {yx a )\ s] = -, Re(s) > |y| 1/a , 

s a — y 


where E a p{x) denotes the generalized Mittag-Leffler (GML) function |22,l23[] defined by 


T,/?U) = ^ 


k =0 


Y(ak + P) 


a > 0, p > 0. 


Then the solution R(x) for m - 1 < p < m (m =1, 2, ...) is written as follows: 

, m -1 


R(x) = -x p E p , p+ i(/) + ^q/ k 1 E p ^ k (x p ). 


(14) 


(15) 


(16) 


(17) 


k =0 
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It should be noticed that at x = 0 the second term of the RHS in Eq. (fl7l) diverges, unless co = 0 for 0 < p < 1, 
and ci = 0 for 0 < p < 2. If p = 1, it is reasonable that Eq. (IT71) reduces to Eq. ©. Therefore, we assume that 
c 0 = ci =0, and we have 


R(x) = -x p E P ' P+ i(x p ). 


(18) 


From Eqs. (fl6l > and (fl8l i. we have 


a —* 


(. x p ) 


k +1 


af^n P (k+ 1) + 1) aU^T(pk+\) 




(x p ) k 




Therefore, Eq. (IT8l > coincides with Eq. © for p = 1. As for the constant Ck, using the following equation, 


oD q x R(x) = -V 

a t—L 


af^n P (k+ 1) + 1) 


h ^ Y -q+p(.k+ 1 ) 1 

,DU p(k+1) = - y r , - n „ „ - -x p -“E p , p - q+l {x p ), 

a ( — n 1 1 1 ' 1 1 ' 


a T(-g + p(k + 1) + 1) a 


which denotes the fractional derivative of R(x) for q > 0 and the fractional integral of R(x) for q < 0, we have 

c k = 0 D p - k -' R(x)\ x=0 = ~x k+ ' E pM 2 (x p )\ x=0 , k = 0, 1. 
a 

Therefore, our assumptions that co = 0 for 0 < p < 1 and co - c\ - 0 for 1 < p < 2 are satisfied. 

Then the solution R(x) of Eq. ( fl8l > and a new form of Planck’s black body radiation law are given respectively 
as follows: 


R(x) = - \E p (x p ) - 1}, x — a(p-Po), 
a/b 


U(x) = 


1 


R(x) EpixP) - 1' 


(19) 


Another possible approach to obtain the fractional blackbody radiation formula is briefly discussed in Appendix B 
Using Eq. d 1 9b . we can analyze NASA COBE data 0], 


3. Various properties of Eqs. ©, Q. ©, and dl9t 

We now investigate various properties of Eqs. ©. Q. ©, and ( | 1 9| i. Concerning Eq. ( | 1 91 ), we have following 
approximate expression with Cb-P o = 0 and x = hv/ksT, and named as (FC II). 


Eq. 03) : f/ (FC U) (T, v, p) = 


C B 

E p (xP) - 1 


-> C B 


1 

- + 

e x - 1 


p — 1 xpi kx k [ij/(k + 1) - In x] 

(e* - 1) 2 ^ r(k + l) 


( 20 ) 


Eq. (l20l) is compared with Eq. 0]). The difference is seen in the second terms of the RHS’s. The behaviors of the 
second terms (named as U( 2 n d)) of the RHS’s in Eqs. ©, ©■ ©, and < 12()b without the coefficients, -//, (q - 1), 
(a - 1), and (p - 1) are presented in Fig. 01 This figure suggests that the Eq. ( IT~9l is almost the same as the 
Bose-Einstein distribution. 

Various analytic corrections for the Stephan-Boltzmann (SB) law without Cg which is expressed as 6/ 2 nd A' 3 
are calculated as follows: 


Bose-Einstein distribution: 


r+oo 3 X 

-Rj 0 = -R- 3 ’ 2 -« 3 ). ( 21 ) 

where /(3) is the Riemann’s 1, function. 

Non-extensive formula II: 


i r 


3 r 2 

x • e x x L 
(e x - l) 2 


dx= (q- 1)-T(6) j . 
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( 22 ) 




















Figure 1: Behavior of the second terms (named as U( 2 nd)) of the RHS’s in Eqs. 0. 0.0. and {J9} without the coefficients factors — p, (q- 1), 
(or - 1), and (p - 1). (a) Eq. 0 (BE), (b) Eq. 0 (NETD II), (c) Eq. 0 (FC I). and (d) Eq. Q9) (FC II) with y3 0 = 0. 


Fractional calculus I: 
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r°° r 3 “ 

fo (e x - l) 2 2 


kx^ij/ik + 1 ) 
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dx = (a 


OO OO , , 


■ 1 ) /1 


T(lr 4 - 11 
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(23) 


Fractional calculus II: 


(P-1) 
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(e x - l) 2 


kx k [if/(k + 1) - In x] 

/ - TTT~. -—- dX 
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OO OO 


<p-» ZZ 
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/ 1 \ 
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i 

- TO + 4) 

In m - - 

+- 

1 I 

\m) 

\k + 3 

k + 2 
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(24) 


Ref. lf26h is utilized in the calculations above. Numerically estimated values of the corrections to the modified 
SB law (U( 2 nd)/C b) are presented in Table [2] It is seen that the roles of // in Eq. Q and (p - 1) in Eq. (l20l) are 
almost the same. 


Table 2: Numerical coefficient factors of fj, (a — 1), (q - 1) and (p — 1) in Eqs. 12 1 1 1241 . 


Bose-Einstein distribution 
Non-extensive formula II 
Fractional calculus I 
(Ertik et al.) 
Fractional calculus II 


-px6£(3) = -px 7.2123414 
(q- l)x 62.22 
( a - l) x 44.41 

(p - 1) x 6.9432884 


4. Analysis of NASA COBE data by Eqs. (Q}, <GK ©• ©. and dlOt 

We are interested in COBE data and now analyze those data in terms of Eqs. (l2Tl> - (l23l i. It is known that 
the COBE data comprise distortion described by a very small dimensionless chemical potential (p) and/or the 
Sunyaev-Zeldovich (SZ) effect IK], and/or a possible effect called the distortion of the space-dimension [27] (see 
the explanation in Introduction.) The first terms of these distributions are the Planck distribution. The second 
terms are different with respect to each other. They correspond to the chemical potential (p) introduced by the 
NASA COBE Collaboration (EH!. 

In particular, our analysis is presented in Fig. [3]and Table[3] Combining the results in Tables[2]and[3] we obtain 
the correction factors for the modified SB law 6'/ Cb for the COBE data. The ratios of the correction factors to 
those of the BE distribution are also shown in the parentheses of Table [4] 
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Figure 2: Behavior of the integrands of Eqs. EMU) without the coefficients factors — p, (q — 1), (a — 1), and (p — 1). (a) Eq. CD (BE), (b) 
Eq. CD (NETD II), (c) Eq. E) (FC I), and (d) Eq. E) (FC II). 



v (cm^ 1 ) 

Figure3: Analysis of the NAS A COBE monopole data in terms of Eq. )20> (FC II): T = 2.7250 K and p- 1 = l.lXlO -5 (x 2 /NDF = 45.0/41). 
Error bars are 400cr. Data are taken from 0 - 


Table 3: Analysis of NASA COBE monopole data in terms of Eqs. {p. CD- E). Q. and ED- 


Eqs. 

T( K) 

/j, (q- 1), (a - 1) or (p - 1) 

^/NDF 

Planck dis. 

2.72502 + 0.00001 

— 

45.1/42 

BE: 0 

2.72501 ± 0.00002 

(-1.1 ±3.2) x 10~ 5 

45.0/41 

NETD II: (0 

2.72502 + 0.00003 

(-0.53 ± 4.98) x lO’ 6 

45.1/41 

FC I: 0 

2.72503 ± 0.00006 

(-0.22 ± 1.38) x 10' 5 

45.1/41 

FC II: 03 

2.72501 + 0.00003 

(1.1 ±3.5)x 10~ 5 

45.0/41 
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Table 4: Corrections to the modified SB law based on analysis of monopole COBE data in terms of Eqs. CD, 0.0.0, and ED- Numerical 
factors in parentheses represent the ratios of the limits to that of the correction for the modified SB law. 


Bose-Einstein dis. 

Non-extensive formula II 

Fractional calculus I 
(Formula proposed by Ertik et al.) 

Fractional calculus II 


-p x 7.2123 = (7.87 ± 23.43) x 1(T 5 
| - p\x 7.2123 < 5.47 x 10' 4 (1.00) 

(q - 1) X 62.22 = (-3.29 ± 31.01) x 1(T 5 
\q — l|x 62.22 < 6.53 x 1(T 4 (1.19) 

(a - 1) x 44.41 = (-9.51 ± 61.01) x 1(T 5 
| ff _ l| x 44.41 < 1.32 x 1(T 3 (2.41) 

(p - 1) x 6.9432 = (7.39 + 24.42) x 1(T 5 
\p - 1| x 6.9432 < 5.62 x 10" 4 (1.03) 


5. Concluding remarks and discussions 


Cl) By applying fractional calculus to the celebrated equation by M. Planck, i.e., Eq. we obtain 
described by the ML function (Eq. (l20l >) which is different from Eq. ([5]) proposed by Ertic et al. 


Appendix A and Appendix B 


a formula 
10. See 


C2) The behavior of Eq. (l20l > is very similar to that of the Bose-Einstein distribution (see Fig. Qjand TableQ}. 

In particular, from analysis of the COBE monopole data in terms of Eqs. Q and (l20l >. we obtain the following 
limits (see TableQ: 


\p\ < 7.6 x 10" 5 (95 %CL), 
\p-l\ < 8.1 x 10~ 5 (95 %CL). 


As is seen in Fig. [4] (residual spectrum (1996)) and Tabled it is difficult to distinguish among the Bose- 
Einstein distribution (f/ (BE) ), Eq. © (j/( NEXTI1 )) based on ^-algebra, and Eq. (l20l > (f/ (FCII) ) based on fractional 
calculus. These are able to describe the distortion of the COBE data. 

The fractional calculus probably reflects the spectral distortion of CMB (Cosmic Microwave Background) of 
the ensemble (the universe) and contains the memory effect at the age of the universe. See also | Appendix B| and 
Appendix C 


Dl) When we analyze the NASA COBE monopole data by means of Eq. (l20t (FCII) including afio = —p, we have 
the following values: T = 2.72502 + 4 x 1(T 5 K, p = (-3.3 ± 5.8) x 1(T 4 , and (p - 1) = (-3.5 ± 6.3) x 1(T 4 . These 
figures depend on the initial values in the CERN MINUIT program. Estimated p and (p - 1) are larger than those 
in Table 0 This fact suggests that we have to choose one freedom between a fit, = -p and (p - 1): In the present 
study a/3 o = 0 i s chosen. The reason of one freedom in FC II is presented in Appendix C Analyzing NASA COBE 


data by Eq. (l20l) including chemical potential p, we find the strong correlation between p and fractional parameter 
(p - 1). This fact implies that we should choose one freedom between them. 


D2) Analyzing the same data by formula of NEXT II, Eq. © with chemical potential p, we find the week corre¬ 
lation between them. See |Appendix C| The magnitude of \p - 1 < 2.3 x 10~ 5 in Table 1C. 51 seems to correspond 
to the SZ effect (y) reported by NASA COBE [32]. (See Ref. (§].) As seen in Fig. 1C.61 b). the situation in FC I, 
i.e., the interrelation between the chemical potential p and fractional parameter (a - 1) is intermediate among three 
cases. 
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Figure 4: Analysis of the residual monopole spectrum, (a) Eq. {2) (BE) with p = -1.1 X 10 5 (x 2 /NDF = 45.9/43). (b) Eq. (NETD II) 
with q - 1 = -0.53 X 1CT 6 (y 2 /NDF = 45.2/43), (c) Eq. 0 (FC I) with a- 1 = -0.22 X 1CT 5 t r/NDF = 46.9/43), and (d) Eq. J20) (FC II) 
with p- 1 = 1.1 X 1CT 5 0 r/NDF = 45.9/43). 


Appendix A. Why do we adopt the fractional derivative in Eq. ©? 


To describe an expansion of the universe at //-era after the Big Bang, Kompaneetz 
propose the following equation for the photon distribution /, 


and Weymann ||29fl 


df 

dt 


kT 6 n e cr e 
m e c 2 c 


dx e 


ox e 


(A.l) 


where T e and T are the electron and radiation temperatures, n e is the electron density, u e is the Thomson cross- 
section and x e = hv/k B T e . The stationary solution of Eq. (IA. 1 b is known as the Bose-Einstein distribution: 


f(x e , //) = 


1 


1 


1 


e x,+p _ ^ 


(A.2) 


where the chemical potential is introduced, reflecting that the number of photons are conserved due to the Compton 
scattering y + e~ y + e~. 

At 380k years after the Big Bang, the CMB (cosmic microwave background) photons are released. The CMB 
photons are described by the Planck distribution, because the number of photons are not conserved. Actually, the 
following processes occur in the universe: The double Compton and Bremsstrahlung scattering y + e~ e + 2y, 
e~ + X —> e~ + X + y, where X denotes an atomic nucleus, usually a proton or Helium-4 nucleus. 

To understand the meaning of the chemical potential /j in Eq. (1A.1I) before the age of recombination in a 
different point of view, we have to take into account the memory effect in Eq. ©: Notice that the variable changes 
to x e —* x — liv/ksT, because the universe has cooled |24:]. 

One of possible methods for taking into account the memory effect of the //-era is an introduction of fractional 
calculus. Following investigations in papers (or books) on geophysics by Caputo and/or Mainardi [3/1 311, we can 
replace the ordinary derivative by the fractional derivative in Eq. ©, i.e., the celebrated equation by M. Planck in 
1900. 


Appendix B. Comparison of the Riemann-Liouville fractional derivative and the Caputo fractional deriva¬ 
tive 

In the present study, we use the Riemann-Liouville fractional derivative to extend Eq. © to fractional order, 
and obtain the fractional blackbody radiation formula, Eq. ( fT~9l >. However, the possible initial condition to Eq. ( 1 1 2| i 
is restricted : As is seen from Eq. (ITTl i. if the initial condition //(()) is not equal to zero, Eq. (IT2l i has no solution. 
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If the Caputo fractional derivative J 30 L [22ll is introduced into Eq. (IT2t instead of the Riemann-Liouville frac¬ 
tional derivative, we can also obtain Eq. ( IT9] > as a solution. Moreover, we can take the initial condition flexibly as 
in the ordinary differential equation. The Caputo fractional derivative of function R(x ) for m = 1, 2, ... is defined 
as 

C 0 D p x R(x) = 1 f (x - T) m - p -'R <m> (T)dT, m - 1 < p < m, (B.l) 

r (m - p) Jo 

where R [m) (T) = (4-) R(t). From Eq. (IB. 1 b . We obtain lim ^D^-R(x) = R lm> (x). The Caputo fractional derivative 

' aT ' p^>m 

is expressed by the use of the Riemann-Liouville fractional derivative and integral as 

{p p x R(x) = (P~ v (pP"‘R(x)), v = m - p > 0. (B.2) 

The difference between the Riemann-Liouville fractional derivative and the Caputo derivative appears distinctly 
in the differentiation of power function x 1 '. If R(x) = x 1 ' with p = 0, 1, m - 1, we obtain that ( ( :D P X x 1 ' = 0 in the 
Caputo fractional derivative. However, x 1 ' = [Tip + 1 )/Y(-p + p + 1) }x~ p+M + 0, in the Riemann-Liouville 
fractional derivative. 

In the definition of the Riemann-Liouville fractional derivative, Eq. m , or the Caputo fractional derivative, 
Eq. ( IB. Il l, integral from x — 0 to x is included. Therefore, the memory effect SHI from the initial stage is taken 
into account in either of fractional derivatives . 


Appendix C. Interrelation between chemical potential p and fractional parameter (p - 1), and the same 
studies for cases of NEXT II and FC I 


To elucidate roles of chemical potential p and fractional parameter (p — 1), we adopt Eq. ( IT9b with afio = p, 

1 


U(x) = 


E p (-(x-p)P) - 1 


(C.l) 


Applying Eq. (1C. 11 to the COBE data, we obtain results in Table 1C.51 To look for an interrelation between p 
and (p — 1), we use a method of Monte Carlo calculus in analyses of COBE data: To understand contents of y 2 


minimum by the CERN MINUET in Table 1C. 51 we adopt an allowed constraint, x 1 ^ (Table 1C.5b + 1.0 B27I1 . 

Moreover, the following set of variables are prepared and 500k generations are performed: 


( T - (T) + 6T x Random number in the [-1, 1] 

| p = (T) + 6p x Random number in the [—1, 1] (C.2) 

( (p - 1) = (p - 1) + 6(p - 1) x Random number in the [—1, 1] 

A number of satisfactory sets (\ 2 < x^ n + L0) is 3021/500k with T = 2.72502 ± 0.00004 K. This ensemble is 
shown in Fg. 1C.51 To investigate an interrelation between (p - 1) and p, we use the method of linear regression. 


Table C.5: Analysis of COBE data by Eqs. (191. Ul and {3} including /?o- 


Eq. 

T (K) 

R 

ip - 1), iq~ 1) °r (a - 1) 

t7ndf 

FC II: O 

2.72502 ± 0.00004 

(-3.3 ± 5.7) x 10~ 4 
\p\ < 1.5 x 10~ 3 (95% CL) 

(-3.5 ± 6.1) x 10~ 4 
\p-l\< 1.6 x 10 -3 (95% CL) 

44.7/40 

NETD II: © 

2.72497 ± 0.00010 

(-3.0 ± 6.3) x 10^ 

|p|< 1.6 xlO- 4 (95% CL) 

(3.5 + 9.7) x 10 -0 

|<7 - 1| < 2.3 x 10~ 5 (95% CL) 

44.9/40 

FC I: © 

2.72492 ± 0.00022 

(-4.4 ± 8.7) x 10^ s 

|p| < 2.2 x IQ" 4 (95% CL) 

(1.5 ± 3.7) x 10~ 3 

\a- 1| < 8.9 x 10~ 5 (95% CL) 

44.8/40 


The following equation is obtained 

(p- 1) = 1.08 p+ 1.04 x 10” 5 (C.3) 

where the correlation coefficient y = 0.998. From Eq. (1C.3l i. it can be said that (p - 1) and p are strongly 
correlated. In other words, (p - 1) and p are not independent each other: The role of the chemical potential p is 
probably replaced by the fractional parameter (p — 1). Indeed this fact is seen in Table|4](in the case of afio — p = 0) 

|p- 1| < 8.1 x 10^ 5 (95% CL) 
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fi (xlO 4 ) n (xlO 5 ) 


Figure C.5: Ensemble of parameters sets (p and (p — 1)) with constraint^ 2 < Xmin (44.66 in Table IC.51 +1.0. a) Number of satisfactory data 
3021/ that of events generations is 500 k. b) Enlarged ensemble with Eq. O- 



Figure C.6: a) Ensemble of parameters sets (p and (q — 1)) with constraint^ 2 < ^ 2 /n (44.86) +1.0. Number of satisfactory data 4269/ that of 
event generation is 200 k. b) Ensemble of parameters sets (ja and (a — 1)) with constraint^ 2 <Xmi n (44.81) +1.0. Number of satisfactory data 
3173/ that of event generation is 500 k. 
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which can be compared with \p - 1| < 1.6 x 10 -3 (95% CL) and \p\ < 1.5 x 1CT 3 (95% CL) in Table [C31 In 
other words, those magnitudes are larger than \p — 1| < 8.1 x 10~ 5 . Thus since we cannot determine p and (p - 1) 
simultaneously in Eq. (Oi l. we choose the fractional parameter (p - 1) between them. 

By making use of the same method for the NEXT II and the FC I, we obtain results in Fig. 1C.61 and the 
following equations: 


(q- 1) = -0.132 p- 5.32 x 1(T 7 (y = -0.859), (C.4) 

(a- 1) = -0.389 a/- 1.72 x 1(T 6 (y = - 0.926). (C.5) 

In the case of NEXT II, the correlation between (q - 1) and p seems to be week. They are probably independent 

quantities, on the contrary to the case of FC II. As seen in Table IC31 and Fig. IC.6l bl. in the case of FC I, the 

situation of (p and (a - 1)) is similar to the case of the NEXT II. 

In conclusion, it should be stressed that those facts above mentioned are based on analyses of COBE data. 
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